Ionic Reactions in Two Dimensions with Disorder 



Jeong-Man Park^^^ and Michael W. Deem^ 
Chemical Engineering Department, University of California, Los Angeles, CA 90095-1592 
^Department of Physics, The Catholic University, Seoul, Korea 

We analyze the dynamics of the ion-dipole pairing reaction in the two-dimensional Coulomb gas 
in the presence of disorder. Sufficiently singular disorder forces the critical temperature of the 
Kosterlitz-Thouless-Berezinskii fixed point to be non-universal. This disorder leads to anomalous 
ion pairing kinetics with a continuously variable decay exponent. Sufficiently strong disorder elim- 
inates the transition altogether. For ions that are chemicaUy reactive, anomalous kinetics with a 
continuously variable decay exponent also occurs in the high-temperature regime. The Coulomb 
interaction inhibits reactant segregation, and so the ionic A'*' + ^ reaction behaves like the 
nonionic A + A ^ reaction. 
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I. INTRODUCTION 

The two-dimensional Coulomb gas has been the sub- 
ject of careful attention since the elucidation of its low- 
temperature phase by Kosterlitz and Thouless [Q and 
Berezinskii |^]. Above a critical value of the dimen- 
sionless temperature, the system approximately obeys 
Debye-Hiickle statistics (as it does for all temperatures 
in three dimensions). Below the transition temperature, 
ions of opposite charge pair to form dipoles. The tem- 
perature at which this metal-insulator transition occurs 
is universal in the absence of disorder, fn the related 
superfluid system, this universality corresponds to the 
universal jump discontinuity in superfluid density (see 
1^ for a review). 

The dynamics of the Coulomb gas under an external 
field has been analyzed by phenomenological extensions 
of the static Kosterlitz- Thouless argument [||^^3j . Differ- 
ent scaling regimes were found, and these are now under- 
stood to correspond to the cases of weak, slowly- varying 
or strong, rapidly- varying external fields Q . While the 
equilibrium properties of the two-dimensional Coulomb 
gas have been established rigorously via field-theoretic 
analysis of the Sine-Gordon Hamiltonian p"5|-p^, there 
has been to date no rigorous, field-theoretic model for the 
ionic dynamics, near the low-temperature critical point 
or otherwise. 

The dynamics of the two-dimensional reaction A"*" -|- 
^ 0, where A"^ and B~ are ions of opposite charge, 
has been studied in the high-temperature limit by scal- 
ing arguments and computer simulation. In the absence 
of Coulomb interaction, the A and B reactants segre- 
gate. This segregation leads to the diffusion-limited de- 
cay law (cA(i)) ~ [no/iSn'^Dt)]'^^'^ Local charge 
neutrality enforced by the Coulomb interaction inhibits 
this segregation of the reactants, allowing for a faster de- 
cay law. The charge density still decays as a power law, 
(cA(i)) ~ at~^. The decay exponent, x, has been ob- 
served in computer simulations to range from 0.79 ±0.04 
^ to 0.85 ± 0.05 Q to unity with logarithmic cor- 



rections [Q. Scaling theories have been proposed that 
lead to values for the decay exponent from 0.85 [^l|j23[ 
to unity [ p4|j25[| . An approximate, self-consistent treat- 
ment of the classical reaction diffusion equations leads 
to the prediction that the decay exponent is unity [^ , 
although logarithmic corrections cannot be excluded due 
to the use of mean-field type equations. 

Studies of single ion diffusion in correlated disor- 
der have shown that for sufficiently long-ranged, disor- 
dered potential fields, anomalous diffusion occurs (see, 
for example, ||ll| , |27| -^ ) . Ionic disorder in two dimen- 
sions creates just such a potential field, w(x), that leads 
to anomalous diffusion. As in previous work, we as- 
sume the potential to be Gaussian, with zero mean and 
correlation function Xwif)- The appropriate form for 
the Fourier transform at long wavelengths is XwO^) = 
J(ixexp(ik • x)x,ji,(x) = l/k"^- This type of disorder 
leads to anomalous diffusion with a continuously variable 
exponent {r'^{t)) ~ bt^-\ where ,5 = 1/[1 + %it / {P'^-f)]. 

In this article, we use the rigorous field-theoretic for- 
mulation of reaction kinetics p5[-p7[ to analyze both the 
ion-pairing reaction near the metal-insulator transition 
and the A^ -I- B~ annihilation reaction at high tem- 
peratures. The master equation description of this re- 
action is described in Section II. The field theory that 
we derive from this description is presented in Section 
III. Two dimensions is the upper critical dimension for 
this system — the dimension below which mean field the- 
ory fails. We derive the renormalization group flows for 
this system in Section IV. We give an asymptotically ex- 
act renormalization group analysis of the long-time dy- 
namics in Section V. For the low temperature phase, we 
find a decay exponent that depends continuously on the 
strength of disorder. Moreover, we find that the critical 
temperature, which is universal in the absence of disor- 
der, depends continuously on the strength of disorder. 
In Section VI we analyze the high-temperature dynam- 
ics of the A"*" -I- B~ — > chemical reaction. We find a 
classical decay in the absence of disorder and anomalous 
kinetics in the presence of disorder. The Coulomb in- 
teraction prevents segregation of the reactants under all 



1 



conditions, and so the dynamics of the ionic A~''+B~ 
reaction is similar to that of the neutral A + A ^ re- 
action. We conclude in section VII with a discussion of 
the experimental implications of our results. 



II. MASTER EQUATION FOR 
LOW-TEMPERATURE ION PAIRING 

To analyze the ion pairing that takes place below the 
transition temperature, we consider the following reac- 
tion 
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where A"*" and B~ are the ions of opposite charge, and 
C is the dipole. We choose initially to have equal den- 
sities of ions (ca(0)) — (cb(0)) = uq and no dipoles. 
The ions are initially distributed at random, with Pois- 
sonian statistics. The long-time decay is not sensitive 
to short-ranged correlations that might be present in 
the initial conditions, such as those resulting from a 
high-temperature quench. The ion-dipole interaction will 
prove to be irrelevant, and so we can ignore the dipole 
orientation. The presence of the dipoles will, however, 
be relevant, and so it is necessary to include the reaction 

®- . . 

By considering the reaction on a lattice, we can write 

a master equation that governs changes in the densities 
of A^, B^, and C. The master equation relates how the 
probability, P, of a given configuration of particles on the 
lattice changes with time: 
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Here mi is the number of A ions on site i, is the number 
of B ions on site i, and li is the number of C particles at 
site i. The summation over i is over all sites on the lattice, 
and the summation over j is over the nearest neighbors of 



site i. The lattice spacing is given by Ar. The diffusive 
transition matrix for hopping from site i to a nearest 
neighbor site j is given by = [1 + (3{u^ — u^)/2] 



and T;^ = [1 + /?(- 



,B „,B 



)/2]. Here u is the sum of an 



external, quenched potential and the Coulomb potential 
created by all of the other ions. Specifically, — Vi + 
Y.k\^i+k - 4o - n,+k]ck and uf = -Vi -f Y.k\^^+k " 
^ko ~ TT^i+k]ck- Here Vi is the external potential at site i, 
and Cfe is the Coulomb interaction c(r) = — Jln(r)/(27r). 
For simplicity we will assume that the ions have the same 
diffusivity, = Db = D. The inverse temperature is 
given by /3= 1/(A:bP). 



III. THE FIELD THEORY 

Using the coherent state representation, we map the 
master equation onto a field theory pSj-pTj. We incor- 
porate a random potential into the field theory via the 
replica trick [ p^j27| . We must also incorporate the ionic 
interaction into the field theory, taking care with self in- 
teraction terms. 

The field theory that we generate is quadratic in the 
fields associated with the dipole density. Integrating out 
these fields, we are left with the action S = So + Si + 
S2 + S3 + S4 + S5 

Sa = J ^'^^JJ dtao,{x,t) [dt - DV"^ + S{t)] a„(x,i) 

d'^x J^' dtba{x, t) [dt - DW^ + S{t)] ha (x, t) 
Si = -no J d'^x [a„ (x, 0) + (x, 0)] 
S'2 = A /" rf'^x / dt[aa{yi,t)baix,t) + aa{yi,t) + ba{ii,t)] 
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Summation is implied over replica indices. The notation 
Jj. stands for / c?'^k/(27r)''. The upper time limit in the 
action is arbitrary as long as it exceeds times for which 
we wish to make calculations. The random, Poissonian 
initial condition is accounted for by the term Si. The 
forward reaction is captured by the term 82- The effec- 
tive potential due to the dipoles is captured by the term 
53. The propagator of the dipoles is given by 



t) 



exphpcfc^ + T-)^], t > 
0,t<0 



(4) 



where Dc is the diffusion coefficient of the dipoles. At 
long times, the ion density will be much smaller than 
the dipole density, and we can replace the instantaneous 
dipole density with the average density. This simplifies 
the effective dipole term to 



^3 — —TlcT 



d'^x / dt 
Jo 

X [aa (X, t)ba (X, t) + da (x, t)+ba{x,t)] , (5) 



with ric — no ~ (cA(i)) ^ no. This modified action is 
identical to that for the reaction 
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This reaction can be recognized as the one addressed by 
the usual Sine-Gordon model of the Coulomb gas, with 
the equilibrium ionic density, y, given at low densities by 

— Uct/X. The fiow equations for these two forms, S3 
and 5*3, are, of course, equivalent. The Coulomb interac- 
tion between the ions is captured by the term 84^. Note 
that the Coulomb coupling should be an effective one, in- 
cluding a finite renormalization due to dipole screening. 
The effective potential due to the randomness is captured 
in the term S5. 

The concentrations, averaged over initial conditions, 
are given by 



(cA(x,t)) = ^im^(a(x, i)) 
(cB(x,i)) = lim {b{x,t)) , 



(7) 



where the average on the right hand side is taken with 
respect to exp{—S). 




The vertices considered in the field theory. We 
have set a* = d + 1 and 6* — 6-1-1. All combinations of 
ip = a,b and — a,b are considered. 



IV. RENORMALIZATION GROUP FLOWS 

We use renormalization group theory to deduce the 
long-time scaling of the ionic concentration. The dia- 
grams that we need to consider are illustrated in Figure 
1. The one-loop flow equations that result are 
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where A is the cutoff in Fourier space, 
exponent is given by 



z = 2- 



(3^j P^-y A2 ncT 
17 ^ I^ZP" AA4 



(8) 

The dynamical 
(9) 



These flow equations are valid to first order in r. At 
this order, they are valid to all orders in pj. Also at 
this order, the flow equation for /?^7 is likely valid to 
all orders in /3^7 The flow equation for A may 

contain contributions from higher orders in /?^7. 



V. MATCHING AND RESULTS 

To compute the long-time value of the ionic concentra- 
tion, we integrate the flow equations up to a matching 
time, to- We match the results of the flow equations to a 
mean field theory that is valid for short times. At these 
short times, we need not worry about renormalization of 
the reaction rates or Coulomb coupling. Furthermore, 
the reaction dynamics occurs in a local region, where the 
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random potential is roughly constant, and so we may 
assume normal diffusive behavior. In other words, we 
can use the standard, classical reaction diffusion equa- 
tions. A self-consistent treatment of these equations has 
recently been presented j2^]. This theory suggests that 
the Coulomb interaction prevents segregation of the re- 
actants. Moreover, the reaction is not limited by local 
transport as long as A < 2pjD. We see that this condi- 
tion is satisfied by the fixed point forward rate, and so 
the concentration is given by 

{cAimj)) = 1/ [l/no{l) + X*t{l)] . (10) 

We find the physical concentration from the relation 

(cA(i)>=e-2'(cA(t(0,0> • (11) 

The result is 

with the fixed point reaction rate given from Eq. (^ as 
X*/D = 2f3J*-l3'^j* = 167r-h/3^7*. Interestingly, we see 
that A* is finite at the Kosterlitz-Thouless fixed point, 
where (tIct)* vanishes. Dipole dissociation, then, is key 
to the physics of the low-temperature fixed point. 

We see that the ions pair according to the classical law 
in the absence of disorder. In the presence of disorder, we 
find anomalous kinetics. The kinetics is anomalous be- 
cause at long-times and low concentrations the reaction 
becomes diffusion limited, and at long times the diffusion 
is anomalous in the type of disorder that we are consid- 
ering. Note that this result for the ion pairing in disorder 
below the transition temperature is identical to that for 
the A -I- A ^ reaction with disorder , except for a 
different value of A* . 

If we interpret these flow equations as relations be- 
tween related equilibrium models, we recognize the stan- 
dard Kosterlitz-Thouless result when disorder is absent, 
with = (ca)^ = ricT/X. Figure 2 shows the flows 
for the case of weak disorder. Interestingly, we can de- 
duce the one-loop critical temperature in the presence of 
disorder with an extension of the elementary Kosterlitz- 
Thouless free energy argument Q|. This argument pre- 
dicts that the ion pairs will unbind when the free energy 
to create two unbound ions, F± — + E-) — r(S'+ + 
5'-), is positive. The Coulomb energy of the ion pair 
is, of course, UG{r) = Jln(r)/(27r). The effective inter- 
action between an ion pair due to the random potential 
is given by Uvcs{r) = -/3~Mn(exp{-/3[w(0) - v{r)]}) = 
— (3j \n{r) / {2tt) . Quenched and annealed statistics are 
identical here for an ion pair seperated by a finte dis- 
tance, r, in a sufficiently large disordered medium, since 
the correlations in the potential for the ion pair are 
short ranged p9| , p0[ |. The entropy of the ion pair is, of 
course, 4A:Bln(r). The ions, therefore, proliferate when 
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FIG. 2. The flow of the equilibrium ion density near the 
low-temperature critical point. Below the critical tempera- 
ture the ion density is driven to zero as dipoles are formed. 
Above the critical temperature the dipoles unpair, and the 
ion density becomes large. For each curve, P^-y = 1 and 
\/D — 16n -f 1 at the point closest to the critical point. 

f3J — /3^7 < Stt. This condition is exactly the one con- 
tained in the flows of Eq. (||) near the low-temperature 
fixed point. This energy-entropy argument is not strictly 
rigorous, since the metal-insulator transition occurs for 
r w L, where L is the system size. In this regime, the 
correlations in the potential for the ion pair are not short- 
ranged, and quenched and annealed statistics are not 
strictly equal. What we have shown is that to one loop 
order these distinct statistics lead to the same behav- 
ior. Unless something unexpected occurs in the regime 
r L, our location of the critical point may be exact to 
all orders. 

The transition temperature, which is universal in the 
absence of disorder, becomes continuously variable in the 
presence of disorder. This is a unique feature of the 
ionic disorder that we are considering. The system un- 
dergoes a transition from insulator to metal either by 
decreasing j3J or by increasing the density of defects, 
p = Figure 3 shows the phase diagram of the 

system at infinitesimally small total (free plus bound) 
ion density. Note that sufficiently strong disorder elimi- 
nates the insulating phase completely. A similar type of 
equilibrium phase diagram has been predicted for two- 
dimensional crystals with random substitutional disorder 
This substitutional disorder is equivalent, in our 
language, to random, quenched dipoles. So we see that 
quenched ions obeying bulk charge neutrality behave in 
the long-wavelength limit in the same way as random, 
quenched dipoles. 

A reentrant metallic phase may occur at low temper- 
atures. This insulating to conducting transition may oc- 
cur because the forces arising from the disorder, which 
tend to separate the ion pairs, are a factor 1/T greater 
than the bare Coulomb forces. Figure 4 shows the reen- 
trant phase diagram predicted by the flow equations 
for p ~ ^Jl/J ~ 0.05 for a range of initial values of 
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VI. HIGH-TEMPERATURE DYNAMICS 
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FIG. 3. The Kosterlitz-Thouless-Berezinskii fixed line in 
tiie presence of disorder (solid) for infinitesimally small total 
ion density. Here p — ^77/ J is roughly the density of defects. 
The system is an insulator below the curve. 



0.030 



0.020 



Conductor 



< 



0.010 



0.000 




1.0 



0.4 0.6 
87t/pj 

FIG. 4. The Kosterlitz-Thouless-Berezinskii fixed line in 
the presence of a fixed amount of disorder, p = y/^/J ~ 0.05. 
Note the reentrant phase at low temperatures for a finite den- 
sity of free ions, y. The curve is strictly valid only in the 
high-temperature regime (solid) . 



y — (ca) = [tIct/A]^/^ and f3J. The temperature at 
which the reentrant phase occurs is roughly proportional 
to ^^7. Since our flow equations are an expansion in 
/9^7, they are not strictly valid in the reentrant regime. 
Thus, the existence of the reentrant phase, while phys- 
ically plausible, cannot be rigorously established with 
our flow equations. A similar reentrant phase diagram 
has been predicted for the equilibrium XY model with 
random Dzyaloshinskii-Moriya interactions, which leads 
again, in our language, to a two-dimensional Coulomb 
gas with random, quenched dipoles Q . 

The ratio p = y^/ J remains constant under rcnormal- 
ization. This means that we can define 7 — 70/e^ and 
J = Jq/sj and one flow equation for e will result. This 
factor e is none other than the dielectric constant! The 
disorder term contains two powers of the dielectric con- 
stant because it is a correlation function of the disorder 
potential. The flow equation for e is not universal 15-1^] 
and should probably include additional (finite) terms. 



We now turn to consider two-dimensional ionic reac- 
tions at high temperatures. That is, we consider the 
reaction 



A+ + B- , 



(13) 



where P is the neutral product of the reaction. In the 
high-temperature regime, the ions pair to an insignificant 
extent. This follows from physical considerations. This 
conclusion also follows from the flow equations that drive 
the ion density to large values. Since the dipole density 
is insignificant, we may ignore the ion-dipole reaction. 
By comparing Eq. ( [l3| ) with Eq. (|^), we see that the 
appropriate action for this reaction is Eq. (^ with the 
replacement A ^ A:, t ^ 0, and ric uq. The flow 
equations for this case are 
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(14) 



In this case, /3J is a constant to all orders. As before 
Jl9| , ^ , it seems likely that /?^7 is a constant to all or- 
ders. The flow equation for k is accurate to flrst order 
only in /3'^7. For our purpose, we will assume that the 
fixed point reaction rate, k* /D = 2(3 J — /3^7, is always 
positive. Note that irrelevant details can renormalize (a 
finite amount) all of the parameters of the model. Dipole 
screening leading to a dielectric constant greater than 
unity is an example of this phenomenon. 

We can again perform the matching. Since at the fixed 
point the reaction step is still rate limiting |26|, we find 
the same classical decay as for the ion-pairing reaction: 



(15) 



Since reactant segregation is suppressed, this result for 
the high-temperature ionic reaction A^ + B^ ^ is iden- 
tical to that for the neutral reaction A + A ^ except 
for a different value of k* 



VII. CONCLUSIONS 

There are many systems well-modeled by the 2-D 
Coulomb gas. A simple physical system might be, for 
example, ions confined to a thin film between two insu- 
lators. Other examples include dislocations or disclina- 
tions in systems such as charge density waves, Abrikosov 
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flux lattices, or Langmuir-Blodgett films. In all cases, 
the defects unbind at higher temperatures, in a form of 
Kosterlitz-Thouless-Berezinskii transition. In the case of 
disclinations, or scalar charges, this transition is exactly 
of the form that we consider, and the system is a per- 
fect instance of the 2-D Coulomb gas model. The type 
of disorder that we consider often comes about in these 
systems via pinning of some of the defects. The density 
of impurities, which are disrupting the low-temperature 
phase, can be controlled via the number of surface defects 
and is given roughly by p = y^/ J- 

For these systems we make the following experimen- 
tal predictions. There should be a continuously- variable 
transition temperature in the presence of long-ranged, 
logarithmic-type disorder. This type of disorder is natu- 
rally induced by impurity phases in these systems. This 
equilibrium behavior has, in fact, been seen in the melt- 
ing of hexatic monolayers and hexatic charge density 
waves 1 45 4^, where disclinations pinned by surface de- 
fects lead to a continuous lowering of the hexatic-liquid 
transition temperature. In other words, these experi- 
ments have shown that the order-disorder transition can 
be driven either by increasing temperature or by in- 
creasing disorder. In terms of Figure 3, these experi- 
ments crossed the transition line by increasing the dis- 
order, i. e. by moving vertically upwards. Ionic reac- 
tions, such as those considered in ||20|-|2^, should decay 
as (cA(i)) ~ l/{2pjDt) at long times in the absence 
of disorder. In the presence of long-ranged, logarithmic- 
type disorder pl] , p7| "p4| , ions at finite density should pair 
in the low-temperature phase according to Eq. (p^). Fi- 
nally, the concentration of ions undergoing a bimolecular 
chemical reaction at high temperature in this same type 
of disorder should decay as Eq. (]l5|). 
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